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We analyze a simple spin-flip process under the presence of two heat reservoirs. While one flip 
process is triggered by a bath at temperature T, the inverse process is activated by a bath at a 
different temperature T' . The situation can be described by using a master equation approach 
in a second quantized Hamiltonian formulation. The stationary solution leads to a generalized 



^ ■ Fermi-Dirac distribution with an effective temperature Tg. Likewise the relaxation time is given in 

terms of Tg. Introducing a spin-representation we perform a Landau expansion for the averaged 
^ I spin (cr) as order parameter and consequently, a free energy functional can be derived. Owing 

' to the two reservoirs the model is invariant with respect to a simultaneous change a ^ — a and 



T T'. This new symmetry generates a third order term in the free energy which gives rise a 
dynamically induced first order transition. 
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I. INTRODUCTION 



Whereas equilibrium statistical mechanics has based on a secure theoretical foundation, this 
is far from being the case in nonequilibrium. The feature of equilibrium phenomena is the 
existence of a probability distribution describing the statistical properties of these systems. 
In general this distribution depends on the interaction among the particles and the tem- 
perature of a single external source called a heat bath. In nonequilibrium the situation is 
different and consequently also the methods in attacking the problems are different. A vari- 
ety of processes are described by Markovian models, where the dynamical process depends 
only on the present configuration of the system. The master equation is one important tool 
for describing different stochastic _grocesses on a complex energy landscape The inputs. 



required for the master equation P| , are a set of states and a set of transition rates between 
those states, for a very recent approach see Often the transition rates are determined in 
according to the principle of detailed balance, in particular in case the system is coupled to 
a single heat bath with a certain but fixed temperature. Thus for stochastic jump processes 
the rates are assumed to follow an Arrhenius ansatz with an activation energy in terms of 
the temperature of the underlying heat bath. Otherwise there is no necessity for having 
only one bath. Therefore we consider here a simple model with two separate heat baths 
with different temperatures. To be specific let us study a annihilation and creation process 
of particles or an equivalent spin-flip process. However both processes should be activated 
by different heat baths. While the spin-flip up-down is triggered by a bath at the temperate 
T, the reversed down-up-fiip is activated by the heat bath at the different temperature T'. 
Apparently both flip rates are likewise determined by different temperatures. One could 
speculate about a generalization by introducing as many heat baths as energy levels exit, 
i.e. each state is related to its own bath and maybe there is an flow between the baths which 
established a typical nonequilibrium situation. 

In the present paper we are interested in a two-level model which can be visualized in terms 
of a spin variable or alternatively by lattice gas variables. The flip process will be orga- 
nized by a coupling to two local baths. An appropriate method to study such a situation 
is given by the master equation approach formulated in terms of second quantized opera- 
tors . In that approach the flip-processes are described by creation and 
annihilation operators, whereas the temperate dependence of the rates are incorporated in 
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the approach by using a Heisenberg-hke picture The approach is generahzed in 

such a manner which enables us to consider two different heat baths. 
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systems with different 
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The analy- 



sis is motivated by searching for some generic features of nonequilibrium steady states. In 
particular, the question arises for a universal behavior under nonequilibrium conditions. As 
an example, a two-temperature, kinetic Ising model is investigated in and extended to 
a diffusive kinetic system in |l6|. The authors found a bicritical point, where two nonequi- 
librium critical lines meet. The analysis is strongly supported by Monte Carlo simulations 
in two dimensions. A similar simulation has been performed studying a two-temperature 
lattice gas model with repulsive interactions . The two dimensional nonequilibrium Ising 
model with competing dynamics induced by two heat baths had been studied in De- 
spite of the two reservoirs the critical exponents belong to the same universality class as 
the corresponding equilibrium model. Alternatively a two temperature lattice gas model 
with repulsive interactions is studied Q, [2^ . Hereby the nonequilibrium transition remains 
continuous unlike in our approach. Another field of interest is the Carnot engine, Carnot 



refrigerator 



19l I27I including a thermally driven ratchet under periodic dichotomous tem- 



perature change |lj], which can be likewise characterized by two reservoirs. General aspects 



22 1 and a rectification 



of a thermodynamic cycle with open flow had been considered in 
of the Clausius inequality is recently discussed in |30]. Totally different physical situations 
occur, when the flow of complex fluids had been analyzed under different heat sources |l7|. 
or in case of a nonlinear oscillator coupled to various heat baths |^|. As pointed out inj^l 
magnetic systems with annealed degrees of freedom are predestined to offer some features of 
a two-temperature systems. An interesting physical explanation for a two reservoir system 
is discussed recently |3| , where the fast and the slow variables of a Hamiltonian system are 
related to different heat baths. In it is demonstrated that the Onsager relations do not 
apply if the two baths are not too close. Apparently transport properties are determined by 
the heat sources. In the occurrence of anomalous heat conductivity in a one-dimensional 
non-Markov process is studied, whereas in 0] a hidden heat transfer is observed, when the 
nonequilibrium steady states are maintained by two heat baths. Very recently in a series 



of papers 



the phase space probability density for steady heat flow is discussed. In that 



case the two reservoirs are mutually connected leading to a flow. 
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As mentioned above we study a spin-flip process under the influence of two heat reservoirs. 
Because this process is on a mesoscopic description related to model A in the classification of 
Hohenberg and Halperin js^ we also analyze the critical dynamics under a two-temperature 
reservoir. 



II. QUANTUM APPROACH TO NONEQUILIBRIUM 

The further analysis is based on a master equation which is written in the form 

dtP{n,t) = CP{n,t). (1) 

Here P{n,t) is the joint probability density that a certain configuration, characterized by 
a state vector n = (ni, ?t,2 . . . jt-at), is realized at time t. In a lattice gas description each 
point is either empty or single occupied leading to rij = 0, 1. Since these numbers can be 
considered as the eigenvalues of the particle number operator and because of the similarity 
of the evolution equation ((T)) to the Schrodinger equation one can introduce a quantum 
formulation of the master equation. This is firstly done by Doi yi for a Bose-like system 
a.d .... by Che. au.ho. fo. spin-ope.... QEffl, .eview L.pa.e gflQ. The 
dynamics of the system is determined completely by the the evolution operator C and the 
commutation relations of the underlying operators. In case of using Pauli-operators the 
restrictions for the occupation numbers to empty and single occupied states is guaranteed, 
see Eq. ©. To transform the basic equation (P) in a second quantized form one has to relate 
the probability distribution P{n, t) to a state vector | F{t) ) in a Fock-space according to 
P(n, t) = (n I F{t)). If the state vectors | n) are a complete set then the last relation implies 
the expansion 

\F(t)) = Y,P{n,t)\n). (2) 

rii 

Under this transformation Eq. (Q) can be rewritten as an equivalent equation in a Fock-space 



d, I F{t)) = L I F{t)) . (3) 

where the operator L is determined in such a manner that its matrix elements correspond to 
C It should be emphasized that the procedure is up to now independent on the realization 
of the basic vectors. As shown by Doi Ijj the average of an arbitrary physical quantity B{n) 
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can be calculated by the average of the corresponding operator B(t) 

{B{t)) = J2Pi^,mn) = {s\B\F{t)) (4) 

rii 

with the state function {s \= "^{n \. Defining the density operator p =\ F{t)){s \ the mean 
value can be even expressed in the conventional manner as 

(S(t))=Tr(p5(t)) . 

The evolution equation for an operator B{t) reads now 

dt{B) = {s\[Bit),L]^\Fit)) (5) 

As the result of the procedure, all the dynamical equations governing the classical problem 
are determined by the structure of the evolution operator L and the commutation rules of 
the operators. In our case the dynamics will be realized by spin-flip processes. 

III. COUPLING TO HEAT BATHS 

Introducing Pauli-operators satisfying the commutation relation 

[d„d]] = 5,,{l-2d4), (6) 
the evolution operator of a flip-process at lattice site i reads fill, compare also Q 

Li = X{dl - did}) + -f{di - d\di) . (7) 
Here the flip-rates A and 7 are parameters which are temperature dependent in case the 
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13l | such a coupling can be 



system is coupled to a heat bath. As demonstrated in 
directly incorporated into the quantum formulation by replacing the operator in Eq. (0) 
through 

L = z/ ^ [(1 - 4) exp{-(3H/2)di exp{(3H/2) + (1 - di) exp{-(3H/2)dl exp{(3H/2)\ . (8) 



The remaining parameter 1/ is determined by the microscopic time scale given by the duration 
of a single spin-flip. The quantity (3 = is the inverse temperature (in units of ks) of the 
heat bath and H describes the static interaction. For a further motivation of this approach 
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and the relation to the Glauber model see also see Now let us generalize the model 
by including two heat baths with different temperatures T and T'. A possible extension of 
Eq. (jSI) is given by 

L = uJ2[i^- 4)e-^^/2^'rf.e^'/2^ + (1 - rf,)e-^^/2^'4e^'/2^] (9) 

i 

The two reservoirs are coupled directly to each lattice point i, therefore the Hamiltonian is 
not the global one but a local energy functional. To illustrate the approach let us discuss 
the simplest case where the Hamiltonian is given by 

H, = {Ei - fi)dldi . (10) 

Here e is a characteristic energy and /i is the chemical potential. Using the algebraic prop- 
erties of the Pauli-operators we get 

L = u^^{l-dl)diexp{{ei- fi)/2T) + {l-di)dlexp{-{e- fi)/"T' . (11) 

i 

Instead of the lattice gas variable rii = d\di with the eigenvalues 0, 1 we can introduce a 
spin variable by CTj = 1 — 2ni. Thus the empty state | ) corresponds to the spin-up state |t) 
and the occupied state | 1 ) is related to the spin-down state ||). Having regard to Eq. (|Tnjl 
we find the non-zero terms of the evolution operator Q 

eM-Hi/2T')diexp{Hi/2T) \ 1) = exp{ei - fi)/2T \ 0) 
exp{~H,/2T')dlexp{H,/2T) \ 0) = exp(-(£, - /i)/2T') | 1) (12) 

The flip process ||) — *|t) is triggered by the heat bath at temperature T whereas the inverse 
process |t) is activated by the bath at T'. Here we have assumed that the activation 

energy e is the same for bath baths. However a generalization to different activation energy 
and consequently different chemical potential is possible. Remark that there is no further 
restrictions for the chemical potential fi Using Eq. © and the algebraic properties of Pauli- 
operators, the evolution equation for the averaged density reads 

u-'dt{ni) = expi-isi - /i)/2T') (1 - m) - expUsi - fi)/2T) {m) (13) 

This equation can be solved exactly and exhibits a stationary solution of the form 



Obviously, the effective temperature is not the mixing temperature of both baths. In a spin 
representation we obtain 

g(£i-M)/2T _ g-(£,-^)/2T' 

In the special case T = T' the stationary solution coincides with the conventional equilibrium 
solution 

(^.). = tanh^ (16) 

If the temperature of one of the heat baths tends to infinity (for instance T' —* oo) the 
stationary solution is 

(a)s = tanh 

When both temperatures T and T' are infinitesimal different from each other T' = T + AT 
the averaged spin is 



AT(£, - /i) 
2r2 sinh(ei - ;u)/T 

The relaxation time tr related to the Eq. p3|) is simply given by 



(a)s = tanh ^ 



(17) 



(z/trei)-^ = exp(^) + exp(-^) (18) 

The relaxation time for T' 7^ T is either enhanced for T' < T or diminished in the opposite 
case. In particular this behavior can be observed for small difference between both baths 
resulting in 



(-t.i)-^ = 2 cosh + ^e~(--)/- 



IV. PHASE TRANSITION 

Now let us study Eq. ()13|) for a magnetic system, where the energy Si at the lattice site 
i depends on the surrounding spin configuration. The chemical potential is in that case 
zero and the interaction is assumed to follow the Ising type. In the simplest mean field 
approximation the energy is given by 

e, = 2[h, + Y,M^3)]. (19) 

i 

where hi is an external field and Jij is the interaction between the z nearest neighbors. 
Firstly the homogeneous case in zero field is discussed. Thus, we have 

e = 2Jz{a) = 2To(a) , (20) 



whereby Tq plays the role of the critical temperature for the conventional case with only one 
heat bath. Inserting Eq. ^ in Eq. ^ we get 
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exp( 



To (a). 



exp 



To±)_ 
T 



exp( 



To(a) . 



(21) 



According to the flip rules defined in Eq. p2|) . this equation is invariant against the simul- 
taneous symmetry transformation a < > —a and T < > T'. Thus, making an expansion 

with respect (a) we find 

iy-'dt{a) = -r{a) + b{af-u{a)\ 
with r = 2Tq 
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(22) 



One can easily check that for T = T' the conventional Landau expansion results. Otherwise 
for two different heat baths with 6 7^ the extended symmetry allows a quadratic term in 
(cr). Thus one concludes that a first order transition is dynamically induced. To illustrate 
the situation in more detail let us derive directly from Eq. ()22j) an analog of the free energy 
for a two bath system: 

_^d{a) OF 



V 



dt 



d{a) 
r 



with 



To + -(o-) 



3^ ' 



(23) 



From here we find another peculiarity of the two bath system, namely the stationary solution 
is different from the minimum of the free energy. The stationary solution follows from 
Eq. (Uni) with Eq. ^ to be 



{a)s = tanh 



(24) 



Apparently this solution is totally different from the minimum of the free energy, which 
follows from Eq. (j22I) by setting dt{a) = 0. 

Now let us consider the case that T' = T + AT. The coefficient, defined in Eq. (j22j) . are for 
temperatures in the vicinity of Tq 



2[t + 



AT, 



AT 



-r- 



u 



h-—] 

3 ^ 4To ^ 



with r 



T-To 



2To^' To' 3^ 4To^ Tq 

When AT ^ T ~ To the coefficient b can be neglected leading to a stationary solution 



(25) 




(To 



AT. 
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In that case a second order phase transition results where the critical temperature is shifted 



to Tn 



AT 



'0 - — • 

In case of an inhomogeneous field hi we calculate the response function defined by 



Xij 



d{a,) 



(26) 



hi=0 



Inserting Eq. (fT!I|) into Eq. (fT!?|l . then the response function fulfills after performing Fourier 
transformation the following equation 

n(T, T') 



with 



m 

n(T, T') 



j(g-)n(T, r) 



n(T, T) 



(27) 



T V 

Here (cr)(t) is replaced by the stationary magnetization which satisfies Eq. fl24j) . In case 
of long wave excitations it is convenient to expand the interaction as 

J(g) = J(0)(1 - eg 2) with J(0) = Jz = Tq . 

Inserting this expression into the solution of Eq. ()27p we get in the limit t ^ oo the stationary 
susceptibility of the form 



with r 



Te - To (1 - as 



In the special case T = T' it results r = [T — Tq{1 — crf )]/(! — erf) ^ T — Tq as expected. 



(28) 



V. CONCLUSIONS 



We have considered spin-flip processes where the rates are conditioned by two different heat 
baths. The situation in mind can be analyzed in a seemingly compact form using the master 
equation in a quantum Hamilton formalism. Using this formalism we find an evolution 
equation for the averaged occupation number, where the stationary solution gives rise to 
a generalized Fermi-Dirac distribution with an effective temperature. This temperature is 
not the mean value of both baths. In terms of equilibrium statistics such a distribution 
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is obtained by coupling separate baths to each energy level. Instead of using a lattice gas 
variable we rewrite the evolution equation in terms of spin variables and end up with a 
Landau-like expansion for the order parameter. Prom here we conclude the existence of a 
free energy functional, where the coefficients depends on both temperatures. Owing to the 
two reservoirs the system allows a new symmetry consisting of the invariance of the evolution 
equation against the change of the spin orientation and simultaneously the interchange of 
the baths. Consequently, a first order phase transition is dynamically induced. In a further 
step both baths should be coupled leading to a temperature fiow. 
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